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ABSTRACT 

Using data from the Sloan Digital Sky Survey Early Data Release and SuperCOS- 
MOS Sky Survey scans of POSS-I plates, we identify a sample of ^2600 subdwarfs 
using reduced proper motion methods and strict selection criteria. This forms one of 
the largest and most reliable samples of candidate subdwarfs known, and enables us to 
determine accurate luminosity functions along many different lines of sight. We derive 
the subdwarf luminosity function with unprecedented accuracy to My < 12.5, finding 
good agreement with recent local estimates but discrepancy with results for the more 
distant spheroid. This provides further evidence that the inner and outer parts of the 
stellar halo cannot be described by a single density distribution. We also find that the 
form of the inner spheroid density profile within heliocentric distances of 2.5 kpc is 
closely matched by a power law with an index of a = —3.15 ± 0.3. 

Key words: stars: low-mass - stars: luminosity function, mass function, subdwarfs 
- Galaxy: halo - Galaxy: structure 



1 INTRODUCTION 

As the most common members of the spheroid, subdwarf 
stars are crucial to understanding the structure and evolu- 
, tion of Galactic halo. Their local scarcity means that they 
are difficult to dete ct, although proper motion studies (eg. 
ICarnev et all fl 9961 have provided several examples of the 
brighter, early-type (F-K) subdwarfs in the solar neighbour- 
hood, yielding important information about their chemical 
composition, density and kinematics. However, only small 
samples exist of both the fainter, later-type subdwarfs be- 
yond the immediate solar neighbourhood. As a result, there 
are significant uncertainties surrounding the characteristics 
of the stellar halo at larger distances and the shape and 
normalisation of the subdwarf luminosity function. 

The advent of large scale, deep surveys with accurate 
photometry such as the SDSS has brought many more of 
these stars within the observational reach of astronomers. 
Determining their number density in the form of their lu- 
minosity function is important to our understanding of the 
Galaxy for a number of reasons. 

Assuming a mass-magnitude relation, the subdwarf lu- 
minosity function can be converted to a mass function, 
the form of which has key implications for both Galac- 
tic and stellar astronomy. On a Galactic level, the sub- 
dwarf luminosity function can be integrated to estimate 
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the number density of subdwarfs in the spheroid. With 
a large sample of subdwarfs from a wide area, the num- 
ber density as a function of line of sight in the Galaxy 
can be used to accurately determine the density law gov- 
erning the spheroid. Additionally, although low mass stars 
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number density can place improved constraints on their 
mass-to-light ratio, and can also be used to attain limits 
of the dark matter contribution of brown dwarfs and white 
dwarfs through extrapolation of the mass function. 



For stellar astronomy, the mass function has important 
implications in understanding the theories of stellar forma- 
tion and evolution. Comparing the mass function of old, 
metal-poor subdwarfs with other types of stars allows de- 
termination of the mass function as a function of epoch and 
metallicity, and hence about the conditions in which the dif- 
ferent populations of stars formed. For example, stars form- 
ing in regions of higher temperature are expected to have 
higher masses. 

As well as providing the stellar mass function, the sub- 
dwarf luminosity function can also enhance understanding 
of globular cluster evolution. Comparing the luminosity and 
mass functions of globular clusters with those of field stars 
of similar ages and metallicities can highlight evolutionary 
processes in the clusters. This method could verify whether 
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the abnormally low mass-to-light ratios of globulars is due 
to evaporation of their low-mass members. 

In this paper we first summarise the status of research 
into the subdwarf luminosity function (Section |2J, before 
describing our data (Section[3J and reduction methods used 
(Sections I4I6H . Finally (Section0 we present our results and 
discuss their consequences. 



2 RECENT RESULTS 

Spheroid members make up only a tiny proportion (~0.2%) 
of the stars in the solar neighbourhood. This factor, coupled 
with their intrinsic faintness, makes subdwarfs hard to de- 
tect, and explains why until very recently samples of known 
subdwarfs have been very small. 

There are two principal methods used to efficiently 
search for subdwarfs. The traditional method is to exploit 
the high heliocentric velocities of spheroid stars by imposing 
a minimum proper motion limit on the sample. This biases 
selection towards spheroid members, increasing the ratio of 
disc to spheroid stars from around 400:1 for a volume-limited 
sample to about 4:1 for the proper motion-limited one. 

The bias this selection introduces can be modelled 
with assumed velocity ellipsoids for the disc and spheroid, 
although there remains a significant sensitivity towards 
the kinematic models used. A further disadvantage of this 
method arises with contamination of the sample from high- 
velocity thick disc stars; however, imposing a strict lower 
tangential velocity cut-off can render this negligible. 

Most studies have employed the method of proper mo- 
tion selection to determine the spheroid luminosity function, 
simply because it remains the most efficient method for ob- 
taining samples of local spheroid stars. Schmidt (1975) used 
it to provide the fir st estimate, with improved determina- 
tions later offered by|Balicall_Jg_Cas^rta no] ll986|) , and many 
in the last decade jPahn et alJ Il995t iGizis fc Reidl 119991 
ICooke fc Reidl200rtlGouldl2003l 1 

The alternative method is to perform deep star counts, 
looking beyond the solar neighbourhood and Galactic disc s 
to probe pure regions of spheroid stars jGould et al.lll998ft . 
This approach has been made possible with the advent 
of sensitive telescopes such as the Hubble Space Telescope 
(HST), enabling accurate star-galaxy separation at large 
distances where the stellar population contains only halo 
members. Whilst samples derived in this way have none of 
the contamination problems of the kinematic studies, there 
are disadvantages. The great distances probed mean that 
trigonometric parallaxes are currently unobtainable, so pho- 
tometric parallaxes are the only indicators of distance and 
hence luminosity. In addition, there is no opportunity to 
obtain follow-up spectra of stars identified in deep studies, 
and such samples can be more prone to bias arising from 
unresolved binaries, although this is countered b y the high 
resolution of HST data in the lGould et alJ Jl99£t) study. 

Despite many investigations into the subdwarf lu- 
minosity function, there remains significant disparity be- 
tween its estimated shape and normalisation . The deep 
star count luminosity function of iGould et al 1 ]l99Sft pre- 
dicts up to five times fewer star s than t hose o f iDahn et alJ 
l|l995l) . lGizis fc Reidl l|l999h and lGouldl l)2003t) for M v > 8, 
and there are even discrepancies among the kinematically- 



selected samples, with the luminosity function of BC86 sub- 
stantially lower than the others. 

There are several possible explanations for the differ- 
ences in the luminosity functions. The first assumption must 
be that the effect is due to systematics: incom pleteness 
in samples (especially iBahcall fc Casertanolll98fitl . inaccu- 
rate photometric parallaxes and differences in kinematic 
models have all been postulated as contributing factors 
{Gould et al.lll998UQizis fc Reidlll999h 

However, the differences between the local, proper 
motion-selected luminosity function s and the deep star 
count result of IGould et al] ]l99St) could partially be a 
manifestat ion of the s pheroi d dens ity profile, a s recog- 
nised by IGould et all dl99gl) and IGizis fc Reidl ll99St) . 
Several studies have indicated that the spheroid has a 
highly flattened centre and a nearly spherical outer com- 
ponent i Sommer-Larsen fc Zhe nll990tlchiba fc Beersll20o"ol : 
IChen. Stouehton fc Smith! l200ll and references therein.) 
This would result in a higher density of spheroid stars for 
the local surveys since they would probe the flattened cen- 
tral ellipsoid while the deeper study would sample regions 
beyond. However, it is unclear whether this structure can 
completely account for the discrepancies between the stud- 
ies, and there are indications that the Galactic stellar popu- 
lations are more complex than those described by the usual 
four discrete components of thin and thick discs, bulge and 
halo. Uncertainties over these associated aspects of Galactic 
structure and discrepancies in the shape and normalisation 
of the subdwarf luminosity function motivate further inves- 
tigation into its true form. 



3 THE DATA 

The sample in this study is drawn from overlapping regions 
of the S loan Digital Sky Surve y (SDSS) Early Data Release 
(EDR) JStoughton et al.l2 002l and scans of P OSS-I E plates 
in the SuperCOSMOS Sky Survey (SSS; iHamblv et ail 
2001al|33-) The areas surveyed are the EDR South and 
North Galactic Stripes (SGS & NGS), outlined in Table □ 
and providing a total coverage of ~394 deg 2 . The SGS over- 
laps 12 SSS fields along the celestial equator, and the NGS 
covers 16 plates (Table Hereafter a 'field' refers to the 
area common to a POSS-I plate and one of the EDR stripes: 
this forms a strip ~2.5°wide along the celestial equator at 

the centre of each plate. 

/ / / / / 

The final u g r i z SDSS filter system is described by 
iFukueita et al] lll996t) . However, at the time of release of 
the EDR th ere were uncertainties in the photometric cali- 
bration (see IStoughton et al 1 I2002T) . so here we refer to the 
preliminary photometry as u* g* r* i* z* . 

For the purposes of measuring proper motions we se- 
lect the POSS-I data for the first epoch (~1950) measures 
and the SDSS as the second epoch (~1998). Whilst proper 
motions are available from the SSS database we make this 
choice so as to maximise the time baseline and hence the 
proper motion accuracy. 

These data form an ideal sample from which to select 
high proper motion spheroid stars: the epoch difference of 
~45 yr enables proper motions to be determined accurately, 
whilst the SDSS photometry provides precise magnitudes 
and colours. 
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Table 1. SDSS EDR regions 



Region RA Dec Area (deg 2 ) 

North Galactic Stripe 145°- 236° 0° 228 
South Galactic Stripe 351°- 56° 0° 166 



Table 2. The survey fields in the EDR SGS & NGS regions 



POSS-I a 8 Plate Area 

Field (J2000.0) 1 b Epoch (xlCT 3 ster) 
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North Galactic Stripe 
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4 METHODS: ASTROMETRY 

The data are first paired, star-galaxy separated, and then 
positional systematics are eliminated by means of an er- 
ror mapping algorithm. Proper motions are derived from 
the two epoch measures, and quasars are used to verify the 
proper motion zero point and to estimate the proper mo- 
tion accuracy. A proper motion cut-off is applied to avoid 
contamination from false detections and to produce a clean 
high proper motion sample. The method of reduced proper 
motion is then employed to select candidate subdwarfs and 
the luminosity function is then derived from this sample. 

4.1 Pairing 

Objects common to both datasets are paired with an algo- 
rithm matching on position, with each object paired to its 
nearest neighbour out to a maximum radius of 10 arcsecs. 
With a mean epoch difference between the POSS-I plates 
and the SDSS data of ~45 yr, this corresponds to an upper 



theoretical proper motion limit of ~220 mas yr~ . Prior to 
pairing the SDSS data are cut at r ^ 20.5, to avoid contam- 
ination from spurious objects on the photographic plates. 
After these cuts there are typically ~50 000 paired objects 
per field. 

4.2 Star-galaxy separation 

The star-galaxy separation uses classifications given in the 
SDSS catalogue. T his has been shown to b e at least 95% ac- 
curate for r* ^ 21 iStoughton et al1l2002l) . adequate for our 
purposes and better than could be obtained from analysis 
of the photographic plates. 

4.3 Position-dependent astrometric errors 

The next stage is to remove position-dependent astrometric 
errors which create well -known 'swirl patterns' on Schmidt 
plates jTaff et aljll992l) of systematic distortions between 
measured plate positions and the expected tangent plane 
coordinates. 

The SSS use s reference stars from the Tycho-2 catalogue 
iHog et alJl2000l) to convert between the measured (x, y) co- 
ordinates and the celestial frame. The reference star cata- 
logue positions are converted to tangent plane (or standard) 
coordinates (g,r?), which are then scale d for Schmidt cubic 
radial distortion jHamblv et al.ll200ld) . Large scale 'swirl 
patterns' are then removed by applying a mean distortion 
map, created by averaging the positional residuals on a grid 
over a large number of plates from the appropriate survey. 
A grid size of 10 arcminute is used since this is adequate 
to map the expected ~30 arcminute scale non-linear distor- 
tions. The mean correction at each location on the plate is 
then obtained by bilinearly interpolating in the grid. These 
corrections are added to the standard coordinates of each 
star, the standard coordinates are fitted to the measured 
positions, and the final conversion to celestial coordinates 
achieved from the reference star plate solution. 

The SSS data therefore has mean large-scale positional 
errors removed, but there will still be plate-to-plate effects 
present as well as systematic errors on smaller scales. These 
are dealt with by an error mapping algorithm first used by 
lEvans fc IrwirJ il995t) (and employed for deriving the SSS 
proper motions), which applies a large scale (10 arcmin or 
lcm) correction to account for plate to plate systematics, 
and then a further one-dimensional algorithm to remove 
small scale (2 arcmin) errors. The procedure operates as 
follows: 

(1) A two-dimensional large-scale error mapper is applied 
by dividing the field into an grid of 10 arcminute squares and 
measuring the mean shift in the x and y coordinates between 
the two datasets from all the stars in each cell. A 3x3 linear 
filter is used to reduce the error function noise, and the 
x and y shift for each object in the field is interpolated 
from the grid point corrections. The SDSS position for each 
object (stars and galaxies) is then shifted with respect to the 
SuperCOSMOS coordinates (the orientation of this shift is 
immaterial) to remove large scale errors. Stars and galaxies 
are not treated separately at this stage so as to facilitate 
the proper motion zero point shift by using the galaxies as 
described in (iv) below. 
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(2) Small scale errors are then removed by a one- 
dimensional error mapper, which calculates x and y errors 
as a function of both x and y by evaluating the mean shift 
for stars in each 2 arcminute strip in both directions over the 
field. Median and then linear filters are applied to smooth 
and reduce the noise in the functions, which are then used to 
again shift the star and galaxy SDSS positions with respect 
to the SuperCOSMOS measures. This process is iterated un- 
til the mean shifts in each of the x and y strips is less than 
a tolerance level (0.008 times the rms x or y shift.) 

(3) The two-dimensional error mapper is then applied 
again to remove any large scale shifts introduced by the 
one-dimensional algorithm. 

(4) At this point the stars have zero mean proper mo- 
tion, whilst the mean galaxy displacement is non-zero, since 
stars dominate the number counts at these magnitudes. The 
galaxies are used to shift the proper motion zero point by de- 
termining a least-squares fit of the galaxy positions to find a 
six-coefficient linear plate model accounting for zero-point, 
scale and orientation. A global translation is then applied 
using this model to reset the galaxies to zero mean proper 
motion, with the SDSS measures transformed by 

xsdss = a + bxsss + cysss 
ysDss = d + exsss + fysss 

where a, b, c, d, e and / are the plate coefficients determined 
from the least-squares fit to the galaxy positions. 

The results of applying this error mapper can be seen in 
Figure n which shows a typical survey field covering a strip 
along the centre of the Schmidt plate that is 6? 25 wide and 
2? 5 high. The first panel shows the field after the SSS mean 
large-scale distortion algorithm has been used, but before 
the error mapping algorithm described above has been ap- 
plied, with the large scale 'swirl pattern' of the photographic 
data still evident. After the error mapper there are only very 
small residual systematic errors in position, and the remain- 
ing random errors have an rms of only <~0.3 arcsec. 



4.4 Photometric-dependent astrometric errors 

Systematic biases dependent on the magnitude of targets 
are particularly prevalent in photographic astronomy, and 
could affect the selection or derived astrometric parameters 
of stars. These systematic magnitude errors cannot be cor- 
rected for in the data, since this could remove real effects 
present, such as brighter stars tending to have higher proper 
motions due to their mean closer proximity to the Sun. How- 
ever, as shown in Figure [5] the influence of these errors is 
likely to be small for the high proper motion sample: their 
size is only <8 mas yr _1 (with a ~45-yr epoch difference) 
for a sample selected with a proper motion cut five times 
larger (see Section r4.6.2L and they have a systematic varia- 
tion with magnitude only at the ~0.1 arcsec level (~2 mas 
yr _1 ). Additionally, this systematic error applies to the error 
difference with magnitude across the whole 15 < r* < 19.5 
range, whereas the subdwarf sample is domina ted by stars 
with r* > 17 llChen et al.ll200ll : iGouldl I2003T) . The more 
pronounced errors at r* < 16 seen in Figure [5] are therefore 
unlikely to have a large effect on the derived proper motions 
in this study, and any magnitude-dependent errors should 
not lead to significant astrometric or selection effects. 
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Figure 1. Before and after the error mapper is applied in field 
0363 in the SGS, which covers a 2?5-high strip along the centre 
of the Schmidt plate. The 'lollipops' show the relative size of the 
binned, smoothed and filtered positional errors at each point in 
the field. Note that the binning and smoothing means that these 
plots show systematic, not random, errors and that the scale of 
the errors in the bottom figure is five times smaller than that in 
the top figure. 



O.i 



u 0.7 

in 

3. 0.6 
o 

ai 0.5 

□ 

■I 0.4 
'<n 
o 

Q_ 

en 0.3 - 
E 

0.2 - 

0.1 


12 




r* magnitude 

Figure 2. The rms of the random positional errors in x and y 
as a function of magnitude for a typical field in the NGS, once 
the error mapping algorithm has been applied. In the magnitude 
range of interest (dashed lines) there is only variation with mag- 
nitude at the ~ 0.1 arcsec level, and the errors are smallest where 
the majority of the subdwarf sample is found (r* ~ 18). These 
considerations indicate that systematic magnitude effects will not 
significantly affect the derived astrometric parameters. 
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Similarly, systematic colour errors are expected to have 
little effect on the astrometry of the sample. Biases aris- 
ing from variations in atmospheric effects are unlikely since 
all observations were taken on or near the meridian, and the 
telescopes used in the survey are at similar latitudes: Apache 
Point Observatory where the SDSS observations were made 
is at 32° N 46', and the Palomar Observatory is at 33° N 21'. 
Discrepancies in positions between the SDSS and SSS data 
due to differential colour refraction will also contribute only 
a small effect, since the SDSS r* and POSS-I E passbands 
are very similar (centred at A ~6400A), and only relative re- 
fraction differences are important when determining proper 
motions. 

4.5 Proper motion accuracy 

4-5. 1 Quasar motions 

The proper motions of all the stars in each field are mea- 
sured after error mapping, and an external check on the 
accuracy of these is obtained from analysis of the quasars 
in each field. The Veron a nd Veron 10th QSO catalogue 
Veron-Cettv fc VeronlEoOlh and the SDSS QSO catalogue 
Schneider et alJ 120021) provide a mean of 135 and 172 
quasars per field in our data for the SGS and NGS respec- 
tively, which are used to estimate the proper motion accu- 
racy. 

The mean spatial discrepancies between the quasars' 
SDSS and SuperCOSMOS coordinates are divided by the 
epoch difference to estimate the rms proper motion error 
for each field, under the assumption that the quasars should 
have zero motion. Due to variations in the quality of the 
plate data this error deviates significantly from field to field, 
but is ~5.0 - 9.0 mas yr -1 for the SGS and ~5.0 - 10.0 mas 
yr -1 for the NGS. However, it should be noted that the true 
proper motion errors of the proper motion sample are liable 
to be less than these values, since the quasars are fainter and 
have poorer centroiding than the majority of the subdwarfs 
in this study. The zero point of the quasar proper motions 
is consistent with zero, as shown in Figure 

4-5.2 Comparison with SSS proper motions 

The accuracy of the derived proper motions can also be as- 
certained by comparis on with those publish ed in the Super- 
COSMOS Sky Survey (lHamblv et al.l2001ch . These are mea- 
sured from the Bj and R plates of the survey, with system- 
atic positional errors removed as in this study (Section 14.31 
and differential colour refraction effects also accounted for. 
The SSS proper motions are currently derived from data at 
two epochs, with a median epoch difference between the Bj 
and R plates of 15 yr, although future releases will soon 
utilise measures from four plates where available. 

To compare the proper motions derived here with the 
SSS results we have analysed SSS UKJ/UKR field 866, 
which overlaps the POSS-1 fields 1440 and 1424 in the NGS. 
The Bj and R plates in field 866 have an epoch difference 
of 17 yr, so with good plate quality this comparison uses 
SSS proper motions that are likely to be of above average 
accuracy. 

The stars in our sample are paired with the SSS data, 
with image quality criteria (such as restrictions on blended 
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Figure 3. Proper motion zero points as measured from quasars 
for the NGS fields (points) and the combined NGS dataset (open 
circle). 

objects, see Section 14.6. H applied as in our analysis. Our 
proper motions derived from the SDSS/SSS data are then 
compared with the SSS proper motions for each star. Fig- 
urcQJshows this comparison for stars with 15.0 < r* < 19.5, 
along with la and 2a deviations from a perfect correlation, 
assuming standard deviations of 10 mas yr -1 and 8 mas 
yr -1 for the SSS and SDSS/SSS proper motions respectively. 
There is good correlation between the measures, with a lin- 
ear correlation coefficient of ~0.85 and ~99% of the stars 
falling with the 2a error bars, indicating negligible system- 
atic effects. Note that there are virtually no stars that have 
discrepantly high SDSS/SSS proper motions compared to 
the SSS measures, although there are some for which the 
SSS proper motion estimate appears to be significantly too 
high. This is likely to be due to the facts that the SSS re- 
sults are derived from a smaller time baseline and from pho- 
tographic material at both epochs, increasing the likelihood 
of objects with spurious proper motions entering the SSS 
sample. The excellent consistency of the SDSS/SSS proper 
motions compared to the SSS measures for stars with high 
SDSS/SSS proper motions demonstrates the good reliability 
of our strict high proper motion sample selection and sug- 
gests negligible contamination from objects with false mo- 
tions. 



4.6 The high proper motion sample 

4-6.1 Image quality criteria 

Prior to the proper motion cut being applied, the data are 
subject to criteria to ensure that only objects with stellar 
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motion selected sample is proportional to the mean popula- 
tion tangential velocity: 



80 120 160 200 

MSDSS/SSS C mos J""" 1 ) 

Figure 4. Comparison of the proper motions derived in this 
study (/^sdss/SSs) with those from the SuperCOSMOS Sky Sur- 
vey (msss) in NGS fields 1424 and 1440. Some ~99% of stars 
fall within the 2a boundaries of a perfect correlation, and the 
linear correlation coefficient is 0.85. The lack of stars along the 
A'SDSS/SSS ^is shows close agreement between the SDSS/SSS 
measures and the corresponding SSS estimates for stars with 
high SDSS/SSS proper motions, demonstrating the good relia- 
bility and negligible contamination of our high proper motion 
sample. 



and high quality images are included in the final sample. 
From the SSS data t he following ch aracteristics of each im- 
age are restricted (see Hamblv et al. 20 01al for more details) : 

• Blend: Objects appearing as blended in the SSS 
database are rejected. 

• Profile Classification Statistic, r\: This quantifies the 
'stellarness' of an object by comparing the residuals of the 
areal profile of an image with that of an average stellar tem- 
plate. Objects with n ^ 4er are rejected. 

• Quality: During processing the SSS data are assigned 
a quality flag, which is affected by circumstances such as 
an image being very large, bright, fragmented or close to a 
bright star or plate boundary. Images with a quality value 
of greater than 127 are rejected. 

• Ellipticity: The ellipticity of an image is calculated from 
the weighted semi-minor and semi-major axes given by the 
SSS processing. Only objects with e ^ 0.25 are included in 
the sample. 



4-6.2 Proper motion selection 

A lower proper motion limit is applied to the sample of 
paired stars for two principal reasons: to amplify the pro- 
portion of spheroid stars selected and to create a 'cleaner' 
proper motion sample. The proper motion selection magni- 
fies the contribution from the higher- velocity spheroid popu- 
lation, since they are effectively sa mpled over larger volumes 
than the lower- velocity disc stars iReidlll99^ : ICooke fc Reidl 
2000). The number of stars of each population in a proper 



N(H > flmin) « po (Vt) 



(1) 



with po the loca l space density of the population llHansonl 
Il985t lReidlll984l) . This therefore amplifies the contribution 
of the higher velocity population above the ratio of the local 
space densities by the amount: 



(2) 



This amplification has a dramatic effect on the likelihood 
of high velocity stars entering the proper motion sample, 
and demonstrates the efficiency of p roper motion se l ection in 
selecting spheroid stars. As seen in ICooke fc Reidl J2000li . a 
spheroid to disc number ratio of Ndisc-N sp heroid = 400:1 for 
a volume-limited sample can be increased to Ndisc-N sp heroid 
= 5:1 for a proper motion sample. 

The second effect of applying a minimum proper motion 
limit to the sample is to include stars with small relative er- 
rors in proper motion, avoiding those with marginal proper 
motions arising from errors in the positions or pairing. This 
results in a 'cleaner' reduced proper motion diagram (Sec- 
tion and a more accurate subdwarf selection. 

The lower proper motion limit can be applied either 
on a global basis, with one limit for all fields, or for each 
field individually, using the proper motion error estimates 
from the quasar positions in each field (Section l4.5ll . Initially 
this latter approach was used, but evidence of systematic 
differences between the proper motion accuracies derived 
for the SGS and NGS led to the adoption of a global proper 
motion minimum for the entire SGS and NGS sample. The 
rms of the quasar 'proper motions' in each stripe was found 
to be a u ~ 6.7 mas yr -1 for the SGS and u M ~ 8.1 mas yr~* 
for the NGS. In order to avoid the low proper motion 'tail' 
in each stripe the more conservative value of the NGS error 
was used for all fields. A 5a u cut on the proper motions was 
applied in defining the lower limit, so that the value used 
for the entire sample was: 



Mmin = 40.5 mas yr 



(3) 



A maximum proper motion limit is also adopted so as to 
ensure sample completeness. This limit is theoretically de- 
termined by the maximum pairing radius of 10 arcsec, cor- 
responding to /i max ~ 220 mas yr -1 . However, the actual 
maximum is somewhat smaller than this value due to sys- 
tematic shifts in position affecting the pairing process; we 
therefore adopt a conservative ^ max = 160 mas yr -1 from 
proper motion number counts (see Section f5. 6. II ) 

Figure |S] shows the relative proportions of high proper 
motion stars with p m in < ^ < (J>max in the SGS and NGS. 
Although there are discrepancies at the bright and faint 
ends, over the magnitude range of the sample selected (15 
< r* < 19) the fraction of stars passing the proper motion 
criteria agree to within a few percent. 
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Figure 5. The numbers of high proper motion stars in the SGS 
(solid line) and NGS (dashed line), expressed as a fraction of the 
total number of paired stars and plotted as a function of magni- 
tude. For the magnitude range of the sample (between the dotted 
lines) the proportion in each stripe agree to within a few per- 
cent, indicating no systematic errors in the derivation of proper 
motions or high proper motion selection. 

5 METHODS: THE SUBDWARF SAMPLE 
5.1 Reduced proper motion 

Although applying a lower proper motion limit produces 
a cleaner sample of fast-moving stars, this will consist not 
only of subdwarfs but also white dwarfs and high velocity 
members of the disc. In order to identify the candidate sub- 
dwarfs, the reduced proper motion of each star is used as a 
discriminator. Defined as 

H — m + 5 log n + 5 = M + 51ogVr - 3.379, (4) 

this parameter separates effectively high- velocity (Vt), low- 
luminosity (M) subdwarfs from the white dwarf and thin 
disc populations. This is seen in the reduced proper motion 
diagram (RPMD) of H r * plotted against (r* —i*): Figures 
HJ and [7| show the RPMDs for all stars with /j, > fimin in the 
SGS and NGS respectively. The fainter and higher tangential 
velocity subdwarfs form a distinct sequence below the thin 
disc main sequence, whilst the white dwarfs populate the 
bluest part of the diagram. 

The colour index (r* — i* ) was chosen because of its ef- 
fectiveness at separating subdwarfs from higher-metallicity 
stars on colour-magnitude diagrams, due to their lower opac- 
ities and higher temperatures at a given mass and luminosity 
resulting in a larger pr oportion of their flux being em itted 
at optical wavelengths JGizislll997| ; lGizis fc Reidll 999). The 
(r* —i*) index is a good temperature/spectral type indica- 



tor for th e M dwarfs which dominate this surv ey (see Sec- 
tion 15.5.11 iLenz et"al ] ll998tlHawlev et aljEoo3) . so on a re- 
duced proper motion or colour-magnitude diagram it enables 
clear distinction between the solar-metallicity stars and the 
hotter (bluer) subdwarfs at each absolute magnitude. The 
(r* — i*) colour is similar to (R — la), proven to be effective 
at separating subdwarfs in this manner, and covers the TiO 
and CaH spectral features important for subdwarf classifica- 
tion jGizisill997llGizis fc Reidlll999t) . We investigated other 
appropriate SDSS indices for use in the RPMD, but (r* — £*) 
proved to be the most effective for subdwarf identification. 

5.2 Reddening and extinction 

With all the fields at high Galactic latitudes (\b\ 37°), 
interstellar extinction and reddening are unlikely to have a 
significant impact on this study. To confirm this we use the 
high-r esolution COBE/DIRBE dust maps of ISchleeel et alJ 
( 1998) to estimate the reddening and extinction at each field 
centre, using the given wavelength coefficients to convert 
Eb-v and Ay to E r *-i* and A r * . We find that the majority 
(23 out of 28) of the SGS and NGS fields have £ r *-i* < 0.04, 
and that the overall mean is £" r *_i* = 0.03. The maximum 
reddening is E r *-i* = 0.08 (for field 0932), and only three 
fields have E r *-i* 0.05, two of which are narrow fields at 
the end of a stripe and hence contain only a small fraction 
of the total sample. In applying these estimates we assume 
that there is no differential reddening across the fields, and 
that all obscuration is in the foreground. This assumption 
is valid for our study, since t he reddening layer scale height 
is ~100 pc iChen et alJl99Sh . and our survey is sensitive to 
very few stars within this distance (see Section f5.5.2l ) 

With photomet ric accuracy of ~3% in the r* and i* 
bands of the EDR dStouehton et al.1 120021) . the size of the 
reddening effect is comparable to that of the oy*-i* « 0.04 
colour errors. Inspection of the reduced proper motion dia- 
grams (Figures[S]and[7J shows that adjustments of this order 
would have little effect on the "tightness" of the population 
locii and subsequent subdwarf selection. The mean extinc- 
tion over all fields is A r — 0.126, which if not corrected for 
would result in derived photometric parallaxes being under- 
estimated by up to ~5%, less than the expected distance 
errors. These considerations suggest that the effects of red- 
dening and extinction on the study are negligible, and so 
corrections are not applied to the data. 

5.3 Tangential velocity cut-off 

Although the RPMD separates the subdwarfs from the thin 
disc stars and white dwarfs, there is some overlap on the dia- 
gram with the locus of thick disc stars. These have relatively 
high velocities, and hence are responsible for the objects ly- 
ing between the old disc and subdwarf sequences. Thick disc 
stars have a much higher local number density than subd- 
warfs, and the inclusion of even a tiny proportion of this 
population in the sample can r esult in a significant overest i- 
mation of the spheroid density feahcall fc CasertandH986l) . 

However, the contamination of the sample by these stars 
can be avoided by the introduction of a cut-off in tangen- 
tial velocity. Figure |H] shows calculated tangential velocity 
distributions of the disc and spheroid populations in the di- 
rection of one of the SGS fields, assuming a velocity ellipsoid 
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Figure 6. The RPMD for stars with fi > ^ min in the SGS. The 
three separate sequences of white dwarfs, subdwarfs and thin disc 
stars are clear. 




1 2 

r*— i* 

Figure 7. The RPMD for stars with u ^ ^ min in the NGS. 



described by a solar motion and Gaus sian dispersion s (Table 
[3J| and following the method given in lMurravl (I1983L p285). 

For our assumed velocity ellipsoids we adopt the disc 
kinematics derived from tw o analyses of local M dwarf s 
jReid. Hawlev fc Gizisl Il995l: lHawlev. Gizis fc Reidl Il99fj) . 
and the spheroid ellipsoid derived from the study of a large 



0.02 




Figure 8. Simulated tangential velocity distributions for the thin 
and thick disc and spheroid (in order of increasing velocity) in the 
direction of field 0363 in the SGS 



Table 3. Ad o pted velocity ellip s oids. The assumed 
disc jR eid et aH |l995t lHawlev et alJ Il996t) and spheroid 
(Chiba & Beers 2000) stellar kinematics used in the tangential 
velocity calculations. All speeds are in kms -1 , and U, V, W 
denote the usual Galactic coordinate triad in the respective 
directions of the Galactic center, direction of rotation and NGP. 
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sample of low metallicity stars bv IChiba fc Beers! {2000). 
The parameters of these ellipsoids are given in Table [3] 

It is clear that a limit of Vt > 200 km s _1 will exclude 
all but a negligible proportion of the thick disc population: 
our calculations suggest that a maximum of just 0.04% of 
the thick disc population will be included in such a sam- 
ple. This cut-off will also cause the low-velocity tail of the 
spheroid population to be excluded from the sample, but 
our calculation results allow us to correspondingly correct 
the derived luminosity functions (see Section [(j.4l ) 

5.4 Selecting candidate subdwarfs 

5.4-1 The colour-magnitude relations 

Equation (QJ suggests that with an assumed colour- 
magnitude relation, lines of constant Vt can be plotted on 
the RPMD. A spheroid-disc separating line of Vt = 200 km 
s _1 can then be applied to the RPMD in order to select can- 
didate subdwarfs, with an upper RPM limit used to exclude 
white dwarfs from the sample. An accurate colour magni- 
tude relation is also required to derive photometric paral- 
laxes for each star in the sample, necessary for computation 
of the luminosity function. 

Obtaining a reliable colour-magnitude relation for sub- 
dwarf stars empirically is problematic, due to the existence 
of relatively few known subdwarfs with accurate trigono- 
metric parallaxes. Theoretical relation s derived from mode l 
atmospheres exist as an alternative jBaraffe et al.llT99^ . 
and these provide excellent matches to cluster sequences 
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and closely follow the majority of field s tar data, espe- 
cially in the infrared fearaffe et alJ Il997l Il998l: IChabrierl 
2003). However, the metal-rich versions ([m/H] > —1.0) 
of these models suffer from discrepancies with the lower 
end of the ob served main sequen c e (V — Ic 2. 0) in op- 
tical colours dBaraffe et al.lll998l : IChabrierl I2003D . Whilst 
lower-metallicity models are expected to be less affected 
in this way, when matched to field stars they still yield 
poor rep roductions of the observed colours a t redder wave- 
lengths fearaffe et alJ Il997l : IChabrierl 120031) . and exhibit 
possible systematic errors in metallicity when used to es- 
timate abundances from colour-magnitude diagrams (iGizisl 
Il997l: iGould et al.lll99l) . These factors could significantly 
affect the ability of the models to provide the reliable opti- 
cal colour-magnitude relations required here, calculating an 
expected absolute magnitude for each subdwarf given an ob- 
served colour and metallicity estimate. Therefore, whilst the 
latest theoretical model atmospheres provide good matches 
to most of the stellar physics and observed parameters, the 
possibility of systematic colour offsets, perhaps wavelength- 
dependent, leads us to choose to employ an empirical colour- 
magnitude relation for this study. 

An obstacle in applying a colour-magnitude relation to 
halo stars is the large range of metallicities present in the 
spheroid, and hence the spread of absolute magnitude for 
a given colour. Without spectra to obtain metallicities we 
deal with this degeneracy in two ways: by using one approx- 
imate colour magnitude relation for selecting subdwarfs on 
the RPMD, and a more accurate estimate which predicts 
metallicity from observed colours for deriving photometric 
parallaxes. 

5.4-2 The colour-magnitude relation for selecting 
subdwarfs 

The colour-magnitude relation used to select subdwarfs 
from the RPMD is obtained in a similar way to that of 
iGizis fc Reidl il999l) . We derive separate relations for low 
and high metallic i ty sub dwarfs, following the classifications 
defined bv IGizisl J 19971) of stars with [Fe/H] « -1.2±0.3 
as subwarfs (sd), and with [Fe/H] ~ — 2.0±0.5 as extreme 
subwarfs (esd). 

To derive a colour-magnitude relatio n we use cal ibrat- 
ing subdw a rfs fro m the compilations of IGizisl il997l) and 
iReid et al.l i200 lT) . Both of these studie s use the spectral 
classification scheme described in IGizisl (Il997ft to identify 
subdwarfs, and give accurate trigonometric parallaxes for 
their samples. These are derived principall y from the USNO 
jMonet et alJ Il992f) and Yale (4th ed; Ivan Altena et alJ 
Il995l) catalogues in the case of Gizis, and from th e Hipparcos 
catalo gue feerrvman et al.lll997t) in the case of IReid et alJ 
feOOll) . 

We consider only subdwarfs with a^/n <0.2 so as to en- 
sure accuracy of the absolute magnitudes, producing a total 
of 35 subdwarfs and extreme subdwarfs from Gizis, and 38 
from Reid et al. We plot these stars on a Mi, (V — Ic) 
Hertzsprung-Russell diagram (Johnson V and Cousins Ic), 
and fit a colour-magnitude relation to the sdM and esdM 
separately. Stellar colour-magnitude relations are poorly 
matched by a linear relation due to an inflection in the sub- 
dwarf main sequence at V-ic ~1.5 (Baraff eet aljll997l) . so 
we match a cubic spline to the parallax data (Figure [5]) 
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Figure 9. Colour-magnitude relations derived from subdwarfs 
with accurate trigonometric parallaxes. Squares denote subdwarfs 
and triang l es ext reme subdwarfs; filled symbols are subdwarfs 
from I Gizisl ll997l) and open symbols are those from IReid et al.l 

l2onih 

There is currently no published SDSS photometry for 
stars with accurate trigonometric parallaxes, so it is nec- 
essary to first derive the colour-magnitude relation in the 
standard Johnson/ Cousins photometric systems, and then 
convert the colours and magnitudes to SDSS r* and r* — i* . 
This is done by means of the relation 

M r - =r*- V + (V-/ C )+M/, (5) 

combined with two-colour relations to convert from SDSS 
(r* - i*) to (r* - V) and (V - ic). 

These two-colour relations should ideally be obtained 
from subdwarfs with accurate photometry and similar 
metallicities to those in our sample. However, of the sub- 
dwar fs in the SDSS Standard Star Catalogue iSmith et alJ 
2002), very few have published Vic photometry and span 
only a small range in colour. We therefore fit linear two- 
colour relations between the Johnson/Cousins and SDSS 
systems for higher-metallicity disc stars instead. Despite 
this, applying these fits to the few subdwarf stars with SDSS 
and Vic magnitudes suggests that the disc metallicity fits 
apply well to the spheroid metallicity stars. The relations 
derived are 

r* — V — — 0.889(r* — i*) — 0.040, (6) 
V - J c = 2.097(r* + 0.429, (7) 

and are shown in Figure Hoi 

5.4-3 Selecting candidate subdwarfs 

With an approximate colour-magnitude relation with which 
to estimate M r * given (r* — i*) for an assumed subdwarf, 
we can make use of Equation Q to plot lines of constant 
V T on the RPMD. 

Figure [TTl stiows the RPMD for the high proper motion 
stars in one of our fields (POSS-I field 0465 in the NGS), 
with the lines corresponding to 200 and 500 km s _1 that 
we use to identify our candidate subdwarfs. These limits 
clearly provide a good match to the observed subdwarf se- 
quence, and are applied on a field-by-field basis to account 
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Figure 10. Two-colour relations from the SDSS to John- 
son/Cousins photometric systems for disc stars in the SDSS stan- 
dard star catalogue. There is little to suggest that these relations 
are strongly dependent on metallicity. 



for the varying kinematics between fields. The lower veloc- 
ity bound is applied to avoid contamination from thick disc 
stars, whereas the upper limit is applied to prevent selec- 
tion of white dwarfs. With knowledge of the expected tan- 
gential velocity distributions for the disc and spheroid stel- 
lar populations (Section 15.31 we can estimate the fraction 
of spheroid and thick disc stars that are expected to fall 
within this velocity range to compensate for selection effects 
("Section lo.41 1. The numbers of high proper motion stars and 
subdwarfs passing all selection criteria (including metallicity 
constraints for subdwarfs) are shown in Table 2] 



5.5 Metallicities and photometric parallaxes 

5.5.1 Estimating metallicities 

The division of the subdwarf colour-magnitude relation into 
two metallicity ranges is clearly a poor approximation, and 
one that is unsuitable for use in obtaining photometric paral- 
laxes of the candidate stars. The spread of nearly two magni- 
tudes in luminosity between metal-rich and metal-poor sub- 
dwarfs at (V — 7c) ~ 2 could lead to significant errors in 
the distance estimates, so a more robust method which ac- 
counts for the absolute magnitude-metallicity degeneracy is 
required. 

The highly accurate five band photometry of the SDSS 
is ideally suited for estimating metallicities from colours. 
The photometric precision is sufficient to clearly separate 
stars of differing metallicities when plotted in metallicity- 
sensitive colours such as (it* — g*) and (g* — r*) (Figure [Toll . 
We use all five photometric bands matched to model atmo- 



Table 4. The numbers of high proper motion stars (AT/,p m ) and 
subdwarfs (N st {) passing all magnitude, proper motion, quality 
and metallicity selection criteria for each field. 
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Figure 11. The RPMD for high proper motion stars in NGS 
field 0465, with the subdwarf (solid lines) and extreme subdwarf 
(dashed lines) selection lines corresponding to Vt = 200 and 500 
kms~ 1 overlaid. 

sphere grids in order to derive a metallicity estimate for each 
subdwarf candidate. 

The model atmos pheres of iLenz et al] |][998) use Ku- 
rucz ATLAS9 models l)Kuruczlll993F to give grids of four 
synthetic SDSS colours for stars with effective temperatures 
(T e ff) from 3500 to 40 000 K and varying metallicity and 
log g. Adopting a value of log g — 4.5 and given the observed 
colours, we use a chi-square statistic and interpolate in the 
grids to estimate the best-fit temperature and metallicity 
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for each subdwarf. Stars with poor metallicity-temperature 
fits are eliminated by rejecting all candidates with a resid- 
ual greater than twenty times the rms value of the residuals 
of the grid fit. This leads to some 3% of candidates being 
rejected, although their highly unusual colours means that 
most of these objects will not be subdwarfs. 

It is clear from inspection of the metallicities derived 
that there is a zero point error in the estimates, with the 
median metallicity for each stripe of [m/H] ~ — 2.4 consid- 
erably _lower_than the [Fe/H]2^j^^5_ecp_ected_for spheroid 
stars jLaird et alJll98StlNissen fc Schusterll99ll) . There are 
a number of likely contributing causes of this: one is that the 
[m/H] = —5.0 model only extends blueward of (g* — r*) ~ 
0.7, and so the metallicities of subdwarfs redder than this 
have to be extrapolated from the [m/H]^ —2.0 models. Even 
subdwarfs with (g* — r*)< 0.7 rely on interpolation between 
the [m/H]= —5.0 and —2.0 models - a very large metallicity 
spread which will inevitably introduce errors. Additionally, 
inaccuracies in the model atmospheres will lead to uncertain- 
ties in the derived metallicities, particularly for redder stars. 
However, the observed offset is largely irrelevant to this ap- 
plication as we are interested in only the relative metallicities 
of the sample. Aside from the zero point error, the derived 
metallicities and temperatures show good consistency, both 
with the expected correlations on a temperature-colour di- 
agram (Figure [T^t and in colour-colour plots (Figure 
Additionally, the metallicity distributions of both the SGS 
and NGS are in close agreement (Figure ITU . indicating no 
significant systematic difference (whether intrinsic or not) 
between the samples. 



5.5.2 Photometric parallaxes 

With a metallicity estimate for each candidate subdwarf, a 
metallicity- dependent colour-magnitude relation is used to 
provide a more reliable indication of each star's intrinsic lu- 
minosity. A median colour-magnitude relation Mi is derived 
by fitting a single cubic spline to all the parallax subdwarfs 
in Section r5.4.2l The absolute magnitude Mj of each star in 
the sample is then estimated by offsetting from this median 
magnitude by an amount proportional to the star's devia- 
tion from the sample median metallicity at that particular 
colour. The estimated absolute magnitude is given by: 



Mi = Mi + d[m/H]w 



where 



dMi 



d[m/H] 



d[m/H] = [m/H]* - [m/H] 



(r*-i*) 



(8) 



(9) 



is the metallicity offset of the star from the sample median 
at that colour, w is a weighting function to allow for colour- 
dependent errors, and the derivative is the variation of ab- 
solute magnitude with metallicity derived from theoretical 
colour-magnitude relations. 

The weight function corrects for the larger colour and 
metallicity errors at redder wavelengths by assigning a 
weight relative to the metallicity spread at each particular 
colour. Hence 
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Figure 12. The correlation between effective temperature and 
r* — i* , with model predictions for [m/H] = —1.0, —2 .0, —5.0 (in 
order of increasing temperature at r* -i» ~ 0.2) from l Lenz et alJ 
1 1998). For the majority of stars (those with r* — i* < 0.4) there 
is close agreement with the models and good internal consistency 
of the metallicity and temperature estimates. The poor fit of the 
models for r* — i* > 0.4 is likely to be mostly due to less accurate 
temperature interpolations over this range, arising from the larger 
colour spread and the fact that the [m/H] = —5.0 model does not 
extend to these redder wavelengths. However, since the majority 
of the sample lies blueward ofr*— i* Ri 0.5 and the metallicities 
are not used in an absolute sense, the effect of these discrepancies 
will be limited. 



where o-*™< e H , is the assumed intrinsic metallicity spread, 
given by the rms metallicity deviation from the median for 
0.1 ^ (r* — i*) ^ 0.25, and oy m / H ^ is the rms metallicity 
spread at a given colour. This weighting therefore assumes 
that the metallicity range for the sample is independent of 
colour. 

The derivative of absolute magnitude with respect to 
metallicity at each colour is derived from the theoretical 
colour-magn itude relations for sta rs of different metallicities 
provided by iBaraffe et alJ II1997I) . The metallicity-absolute 
magnitude relation varies considerably with colour (Fig- 
ure but is approximately linear. The slope of this rela- 
tion is therefore a constant and is estimated from the mod- 
els for each star's colou r to provide the last term in Equa- 
tion JHJ. The models of Ba raffe et alJ (|l997) and the paral- 
lax subdwarfs are given in the {Mi, V — Ic) colour plane; 
conversion to SDSS colours is achieved using Equation @. 

It is stressed that this modification of the colour mag- 
nitude relation based upon estimated metallicity is only an 
approximation. It is only intended to adjust the predicted 
magnitude of each star according to its likely metallicity in 
the right direction and by roughly the right amount. In light 
of this the adjustment of absolute magnitude from the me- 
dian is restricted to a maximum of ±1.5 mag, to prevent 
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Figure 15. The correla tion of absolu t e mag nitude with metal- 
licity from the models of lBaraffe et alj Il997l) . for (V — Iq) from 
1.0 (top) to 2.8 (bottom) in steps of 0.2. 



Figure 13. Colour-colour plots for subdwarfs in the SGS (top) 
and NGS (bottom), with model atmosphere predictions for stars 
of [m/H] = +1.0,0.0,-1.0,-2.0 and -5.0 (from bottom to top) 
overlaid. There is no evidence for any systematic metallicity dif- 
ference between the samples. 
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Figure 14. The estimated metallicity distributions for the SGS 
(solid line) and NGS (dashed line). Although the metallicities 
are clearly not accurate in an absolute sense from comparison 
with the [Fe/H] ~ —1.5 expected for spheroid stars, there is good 
agreement between the samples and no systematic metallicity dif- 
ference between the two. The median metallicities of —2.5 for the 
SGS and —2.3 for the NGS agree well within the expected errors. 



any stars being assigned inappropriate values. However, it 
can be seen from the sample estimated colour-magnitude 
diagram (Figure [Tol that this method provides a far more 
realistic and accurate distribution of absolute magnitudes 
than a simple division of the colour-magnitude relation into 
two metallicity ranges. 

With an estimated absolute magnitude for each star a 
distance can then be derived, assuming no significant red- 
dening effects due to the high latitude of the fields (\b\ ^ 
37°). The distances sampled are in the range ~260 pc to 
~2.8 kpc, and the distributions for the SGS and NGS are 
shown in Figure 1171 That proportionally fewer stars with 
distances over 2 kpc are found in the SGS compared to the 
NGS is an effect expected in a non-uniform spheroid den- 
sity distribution. Whereas the NGS is directed principally 
towards the inner quadrants of the Galaxy, the SGS lies to- 
wards the outer quadrants where the space density of stars 
will be lower with a non-uniform distribution. 



5.6 Tests for sample systematic errors 

The luminosity functions derived from the two independent 
SGS and NGS samples not only provide an excellent probe of 
the spheroid density distribution, but also can be compared 
to ensure that our results are free from significant systematic 
errors - a source of concern with earlier studies (Section 

Prior to comparing the results from each stripe, how- 
ever, we must ensure that the two samples are self-consistent 
and free from systematic effects. Inconsistency between the 
samples could arise from either a systematic disparity in the 
selection of each sample, or from an intrinsic difference in 
the constituent stars of each. 
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Figure 16. Estimated absolute magnitudes for all subdwarfs in 
the SGS and NGS, derived from a median colour-magnitude rela- 
tion fitted to calibrating subdwarfs (solid line) and offset accord- 
ing to the metallicity estimate of each star. 



5.6.1 Sampling systematic errors 

Systematic errors arising from the sample selection are the 
most likely, with a whole range of possible causes such as a 
difference in astrometric or photometric accuracy between 
the stripes. However, we can test for any significant differ- 
ence of this type by investigating the completeness of each 
subdwarf sample. We do this in a somewhat crude manner 
by comparing cumulative proper motion number counts: as- 
suming a uniform stellar density and that proper motion 
and distance are inversely proportional (/x cx cP 1 ), a plot 
of log cumulative number count (SJV) against log [i should 
have a gradient of —3: 



E7V cx V <x d 3 cx fi~ 3 . 



(11) 



Figure 1181 shows the histogram of cumulative proper 
motion counts for all paired stars in the NGS and SGS after 
the error mapper has been applied. We fit a straight line to 
the points between the global proper motion limits for each 
stripe: fjmin = 40.5 mas yr _1 and (i ma i = 160 mas yr _1 . 
Both give a close fit to the expected gradient: — 2.92± 0.06 
for the SGS and -2.99± 0.06 for the NGS. Similar plots 
for just the candidate subdwarfs in each stripe are shown 
in Figure ED where gradients of —3.10 are found for each 
stripe. 

There is some evidence for incompleteness in the non- 
linearity of the histograms, especially for the subdwarf sam- 
ples where proper motion errors lead to incompleteness to- 
wards the proper motion limits. This is particularly relevant 
to the higher proper motions where the errors will be more 
dominant, and this explains the rapid tailing off of the distri- 
butions for log n > 2.1 caused by the imposition of the upper 
proper motion limit. However, the number counts for the 
population from which the subdwarfs are selected (Figure 
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Figure 17. The distributions of estimated distances for the SGS 
(top) and NGS (bottom). 



1181 demonstrate that it is complete for 40.5 mas yr _1 ^C /i ^ 
160 mas yr~ , and the good match of all of the histograms to 
the expected gradient of —3 indicate that there are no signif- 
icant signs of incompleteness affecting the samples. It must 
be also stressed that these are only approximate tests for 
completeness, and will have some unreliability introduced 
by the assumption of uniform stellar density, which is par- 
ticularly invalid for stars at large distances which tend to 
have lower proper motions. 



5.6.2 Radial metallicity gradient 

An intrinsic difference between the samples that could have 
a significant effect on the results is a difference in median 
metallicity between the SGS and NGS. Although the radial 
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Gradient (1.61.2.20): -2.923 ± 0.060 
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Gradient (1.61,2.20): -2.988 ± 0.062 
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Figure 18. Cumulative proper motion number counts for all 
stars in the SGS (top) and NGS (bottom). The straight line fits 
between the proper motion limits in each stripe are shown (dashed 
line), along with the expected line of gradient of —3.0 (dotted 
line). 
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Figure 19. Cumulative proper motion number counts for can- 
didate subdwarfs in the SGS (lower line) and NGS (upper line). 
The straight line fits between 1.6 ^ log/x ^ 2.1 (dashed lines) 
both have gradients of —3.1 



the luminosity function would be underestimated. The de- 
rived luminosity function of each sample would therefore be 
expected to differ. This effect would be less pronounced with 
the photometric parallax relation defined in Section 15.5.21 
that partially accounts for metallicity variations, but is an 
issue that should still be addressed. 

However, our data display little evidence of metallicity 
gradient. The two-colour diagrams for the subdwarf sample 
in each stripe in Figure lTHl show no systematic difference, and 
the estimated metallicity distribution and median metallic- 
ity of each (Figure tTH are very similar. It is therefore highly 
unlikely that there is any strong systematic metallicity dif- 
ference between the samples. 

These checks of sample consistency indicate that there 
are no large systematic differences between the samples aris- 
ing from sampling errors or intrinsic variations. A further 
and final test of the sample completenesses can be ascer- 
tained following the luminosity function calculations from 
the value of (V/V max ) (Section 01. 

6 METHODS: THE LUMINOSITY FUNCTION 

6.1 The generalised Vmax method 

With a final sample size of 918 candidate subdwarfs in the 
SGS and 1696 in the NGS, the luminosity function can be 
accurately derived. We achieve this using a modification of 
the Vmax method of ISchmidtl il96ct) : each star contributes 
l/Vmax to the luminosity function, where V ma x is the maxi- 
mum volume that the star could have been detected in, given 
the proper motion and magnitude limits of the survey. This 
technique therefore implicitly corrects for any bias arising 
from the proper motion selection. 

Schmidt's original 1/Vmax method assumes that the 
sample is selected from a uniformly-distributed popula- 
tion. This certainly is not the case for our subdwarf sam- 
ple; indeed, we intend to use it to determine the vari- 
ation in spheroid densit y. We therefore adopt the 'gen- 
eralised Vmax' 1 method llStobie. Ishida. fc Peacockl Il989l : 
iTinnev. Reid. fc Mouldl Il993h . which extends Schmidt's 
method to non-uniformly distributed samples. Given that 
each star in the survey can be detected to minimum and 
maximum distances d m i n ,d ma x in a field of solid angle f2, 
the luminosity function is defined as: 

^sample = ^ 7F—> ( 12 ) 
V max 

vLax = Sl [ ""Vn(£) dfl) (13) 

Jd min n © 

where s is the heliocentric distance and n(s) is the adopted 
spheroid density law (Section 16. M^ , 



metallicity gradient of the inner spheroid i s thought to be 
small (A[Fe/H]/Ar w -0.012 dex kpc" 1 . iBekki fc Chibal 
l200ll) . with both samples in opposite radial Galactic direc- 
tions this could nevertheless introduce a source of error. The 
assumption of a single colour-magnitude relation for both 
of the samples would mean that stars in the sample with 
the higher metallicity would have their absolute magnitudes 
overestimated, and hence their distances and contribution to 



6.2 Proper motion and magnitude limits 

The distance limits depend on the minimum and maximum 
proper motion and magnitudes in each field. The upper and 
lower proper motion limits employed are fimin = 40.5 mas 
yr _1 and fj, max = 160 mas yr _1 (Section l4.6.2t . 

Ascertaining the magnitude limits for each field is less 
straightforward as the data come from two quite different 
sources, and there is significant variation in the depth of 



The Subdwarf Luminosity Function 15 




Figure 20. The log number - magnitude histogram for all un- 
paired (top) and paired (bottom) stars in the NGS, with the 
approximate upper and lower r* magnitude limits for the NGS 
fields shown. The actual magnitude limits adopted for each field 
are ascertained from where the log(A r (r*)) curve in the field's 
histogram departs from a steady increase. 



the plate material. However, the fact that the datasets are 
paired and that the SDSS data probe much fainter (r* ~ 23) 
than the POSS-I plates means that only the r* magnitudes 
need be considered. 

The POSS-I plates reach incomplete levels at R ~ 20, 
at which magnitude the SDSS data are certainly complete. 
Stars need to appear in both datasets to pass the pairing 
criteria, so the faint magnitude limit can be defined solely 
in terms of the much more accurate SDSS r* magnitude. A 
(log number, r* magnitude) histogram is plotted for all of 
the paired stars in each field, and an upper limit is conser- 
vatively defined from the magnitude at which log (JV(r*)) 
departs from a steady increase (Figure 1201 Table |K|) The 
brighter limit is less crucial due to the smaller likelihood of 
subdwarfs having such magnitudes, but a limit of r* ~ 15 is 
applied to encompass the range that is likely to be complete 
as indicated by the histogram. 

With proper motion limits /Umin, Mmax and magnitude 
limits r^ in , r^ax for each field, the limiting distances for 
a subdwarf of distance d and magnitude r* are therefore 
defined as 



rfmin = d max 



d min 



10' 



,0.2(r 



/^rnax 



10' 



(14) 



(15) 



6.3 The spheroid density law 

Recent studies of the structure of the Galactic spheroid have 
indicated that its density profile is flattened and follows a 



puwei law: an axial lanu oi ^c/a 

ilGould et alJll998HSluis & Arnold! 


i«s u.u aim fJ\ i } ~ r 
19981: lYanny et al ll200ot 


Ilvezic et alj2000HChen et al.l200li 


Sicsel et al.l2002HGould 


2003). We therefore adopt a spheroid density law of the form 



rl + Rl 



n© 



+ R 2 + z 2 /q 2 



Table 5. Magnitude limits and spheroid and thick disc discovery 
fractions for SGS (top) and NGS (bottom) fields. The minimum 
magnitude limit for all fields is r* = 15.0 and the proper motion 
limits are /i m i n = 40.5 mas yr — 1 and fimax = 160 mas yr — 1 . 



Field 


r* 

i'l'i o. ,r 


Xh(%) 


Xtd{%) 


0932 


18.7 


61.15 


0.01 


0363 


19.5 


62.56 


0.01 


1453 


19.0 


63.57 


0.01 


1283 


19.2 


64.16 


0.02 


0852 


19.1 


64.39 


0.02 


0362 


19.2 


64.27 


0.02 


1259 


19.1 


63.85 


0.03 


1196 


19.1 


63.13 


0.03 


0591 


19.3 


62.13 


0.04 


0319 


19.2 


60.85 


0.04 


0431 


19.0 


59.31 


0.04 


0834 


18.9 


57.50 


0.04 


0151 


19.0 


61.15 


0.01 


1402 


19.2 


62.56 


0.01 


1613 


18.8 


63.57 


0.01 


1440 


19.1 


64.16 


0.02 


1424 


19.1 


64.39 


0.02 


0465 


19.0 


64.27 


0.02 


1595 


18.7 


63.87 


0.03 


1578 


18.9 


63.13 


0.03 


1405 


19.0 


62.12 


0.04 


1401 


19.2 


60.85 


0.04 


0471 


19.2 


59.31 


0.04 


1400 


19.0 


57.50 


0.04 


1397 


18.7 


55.42 


0.03 


0467 


19.3 


53.26 


0.03 


0470 


19.1 


50.94 


0.03 


1318 


18.8 


48.62 


0.03 



where n© is the local spheroid density, a is the power law 
index and q — (c/a) is the axial ratio, and we assume Rq 
= 8.0 kpc and r c = 1.0 kpc. Converting the Galactocentric 
cylindrical coordinates R and z in terms of heliocentric co- 
ordinates (s,l,b), where s is the distance from the Sun and 
I and b are Galactic latitude and longitude: 



+ Rl 



n© 



+ Rl + s 2 



2sR cos I cos b + s 2 sin 2 b^jf- 



.(17) 



(16) 



The two density law parameters a and q are allowed 
to vary, and luminosity functions are derived for a range of 
values to find the (a, q) law which best matches the observa- 
tions (see Section fo-SI ) A significant benefit of this approach 
is that the density law variables a and q can be fitted with- 
out having to make an assumption about the local spheroid 
density tiq. 

6.4 Discovery fraction 

Given distance limits and an assumed (a, q) density distribu- 
tion, the sample luminosity function can then be calculated. 
However, the tangential velocity cut-offs cause the sample 
to exclude a given fraction of spheroid stars, so this scaling 
needs to be taken into account in order to derive the true 
spheroid luminosity function from the sample. 

The fraction \h of spheroid stars expected to have 200 
km s _1 ^ Vt ^ 500 km s _1 is calculated for each field using 
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tangential velocity calculations as described in Section [5.31 
Using our adopted velocity ellipsoids, \h varies between 0.58 
and 0.64 in the SGS, and 0.49 and 0.64 in the NGS (Table0, 
and scales the spheroid luminosity function as 



*■ ^sample 

Xh h 



(18) 



At this stage any possible contamination by thick disc 
stars can be considered. Assuming they are also included in 
the sample with Vr ^ 200km s , then the derived lumi- 
nosity function will comprise a total for spheroid and thick 
disc members. The sample spheroid luminosity function can 
then be calculated from the total (spheroid plus thick disc) 
by 



4> 



sample 



sample 
b+td ' 



(19) 



where Xh is the fraction of spheroid stars in the sample. This 
is given by 



X h = 



[Xtd /Xh)(ntd/n h ) + l' 



(20) 



where \td and Xh are the fractions of thick disc and 
spheroid stars with 200km s _1 ?C Vr ^ 500 km s _1 
and ntd and rih are the local number densities of thick 
disc and spheroid stars. The discovery fractions xtd and 
Xh are known from the calculations described in Sec- 
tion 15.31 whilst the relative normalisation of thick disc to 
spheroid stars is taken from independent studi es. Assum- 
ing a thick to thin disc density ratio of 1:10 ||Rej clet_alJ 
119951: IChen. Stoughton fe Smithl 1200 it ISieeel et alJ 120021) 
and a combined disc to spheroid normalisation of 400:1 
JChabrier fc Merall997llHormberg fc Flvnrl2000l) . we adopt 
a thick disc to spheroid ratio of n t d:rih = 40:1. With this 
consideration of thick disc contamination, the true spheroid 
luminosity function is therefore derived from 



Xh 



sample 
td * 



(21) 



However, with a tangential velocity cut-off of Vr ^ 200 
km s _1 , this scaling for thick disc contamination has very 
little effect on the luminosity function. This normalisation 
gives a scaling factor of 0.97 ^ Ah ^ 1.00, so at worst the 
thick disc contamination has just a three percent effect on 
the normalisation of the luminosity function. 



6.5 Luminosity function errors 

In estimating the errors in the lumin osity function we adopt 
the assumption of Poissonian errors l)Feltenlll97ri : 



2 



\ " 



p 1c /L^i V 2 

V max 

Allowing for scaling: 



2 / Aft \ 2 



6.6 Combining fields 



(22) 



(23) 



Whilst the luminosity functions and densities for each field 
are derived separately to investigate the spheroid density 
profile, it is desirable to calculate a combined luminosity 



function for the fields in each stripe. Unfortunately the ne- 
cessity of scaling to account for the spheroid discovery frac- 
tion and thick disc contamination on a field by field basis 
means that a total 1/Vmax luminosity function cannot be 
calculated for the whole sample. However, a combined lu- 
minosity function for all of the fields in each stripe can be 
derived by combining the luminosity functions for each field 
with a simple weighted mean. The mean density and error 
for each magnitude bin are then given by 



log $ = 



i 



(24) 
(25) 



where the summations are over the fields in each stripe. 

6.7 Transforming $(M r «) to $(M V ) 

To facilitate comparison with published luminosity func- 
tions, it is necessary to convert the derived luminosity func- 
tion from SDSS r* absolute magnitude to Johnson V . This 
is achieved using the transformation 

dM r * 



HMv) = HM r *, dMv . 
The derivative is evaluated from 



Ah 



M r * 

M r * 



(r* - V) 
[(-0.889(r* 



0.040] , 



(26) 



(27) 
(28) 



from Equation (JSJ. By varying (r* — i*), M r * can be plot- 
ted against Mv, a spline fitted and the derivative evaluated 
(Figure|2TJ- Hence dM r * /dM v can be calculated for a given 
M r * , and repeating this for each M r » bin in the luminosity 
function completes the transformation. 



6.8 Constraining the spheroid density profile 

With a wide range of lines of sight, the SGS and NGS sam- 
ples are ideally suited to constraining the density distribu- 
tion of the spheroid. Indeed, a non-uniform density law has 
to be assumed in order to compare the luminosity functions 
of the different samples. Figure l2~2l shows the number densi- 
ties (luminosity function integrated over 5 ^ M r » ^ 10) for 
each field plotted against Galactic longitude when derived 
using a conventional V m ax method under the assumption of 
uniform space density. It is clear from comparison with the 
models that at the distances sampled the varying density of 
the spheroid has a strong effect and so has to be taken into 
account when deriving the luminosity function. 

In order to constrain the density distribution a power 
law of the form given in Equation 1171 is assumed, and sub- 
dwarf number densities are calculated for each field in the 
SGS and NGS for a range of power law indices and axial 
ratios. The power law index a is allowed to vary from —2.0 
to —4.4 in steps of 0.05, and the axial ratio from 0.2 to 1.0 
in intervals of 0.05. 

The generalised V ma x method (Section l6.H should en- 
sure that the derived number densities are constant and in- 
dependent of line of sight, so how well a (a, q) power law 
model fits the data is ascertained by measuring how close 
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1.05 



Figure 21. The relations between M r * and My (solid line) and 
d J^f^ (dot-dashed line) derived from the two-colour relations. 
These are used to convert the luminosity function derived in M r * 
to My to facilitate comparison with published results. 



to a uniform value the field number densities are under the 
model. This goodness of fit is defined by a modified chi- 
square statistic 



2 

X = 



EE 1 ? 



= 1 3 = 



2 +<7 2 
» 1 J 



(29) 



where n, and a; are the number density and its standard 
deviation for each of nf fields. The best-fit (a, q) density 
model is therefore the model with the minimum value of 
X 2 , and the likelihood of each model can be evaluated by 
comparing its statistic value with the minimum. 



7 RESULTS AND DISCUSSION 
7.1 Spheroid density profile 

Figure 1231 shows the results of fitting spheroid density mod- 
els (a, q) to the combined SGS and NGS data, with con- 
tours of equal \ 2 plotted. Whilst the axial ratio q can- 
not be constrained with these data, the power law in- 
dex a can. Although there is a slight degeneracy with 
q, the best-fit power law index is a = — 3.15± 0.30 for 
the range 0.55 ^ q ^ 0.85, which is t he value of the 
axial ratio derived from recent studies <|Sluis fc A rnold 
1998tlGould. Flvnn fc Bahcalllll99cl iRobin. Revle. fc Crezel 
200ol: IChen. Stouehton fe Smith) 120011: ISieeel et alJ 12002^ . 



This power law index is largely in agreement with the re- 
cent results of iGould et al.1 Jl998t) (a = - 3.13 ± 0.23); 
ISluis fc Arnoldl (Il998l) 7a = -3.2 ± 0.3): [Yannv et alJ 
fcOOOl) (a = -3.2 ± 0.3) ; llvezic et alJ foOOfj) (a = -3): 
Chen. Stouehton fc SmitU ll200ll) (a = -2.55 ± 0.3); 
Sieeel et alJ ll2002l) (a -2.75) and iGouldl i2003h (a = 




-3.1 ± 1.0). This r distribution differs significantly from 



-30 30 

Galactic Latitude 

Figure 22. The subdwarf number densities (points) in each field 
derived from a conventional V m ax method assuming uniform pop- 
ulation density. The lines show the expected relative number den- 
sities at a distance of lkpc from a spheroid density power law with 
(a, q) = (—3.0, 0.55) . Is is clear that the uniform density assump- 
tion is invalid for this sample and so a non-uniform distribution 
must be assumed in the calculation of the luminosity functions. 
Accordingly, these data can be used to constrain the form of the 
density profile. 



the r~ 2 power law of the dark matter halo, providing fur- 
ther evidence for the theory that faint halo stars consti- 
tute iat_jnost L _a_negliribl£_j^ 

ter jBahcall et alJl99llFlvnn et alll99dlc"habrier fc Meral 
ll997tlFields et alJll998l) . 

These results are only able to constrain q > 0.3. Com- 
pariso n of the limits of q ~ 0.55 ± 0.06 set bv lChen et alJ 
( 2001) with the SDSS EDR data and inspection of Equa- 
tion 117H shows that this is due to the much smaller distances 
probed by this study. Small values of s in Equation 1171 
cause the q term to be negligible compared to the r 2 + R 2 
terms, and so variations of I and b are insensitive to q. The 
low distances are a result of the relatively bright lower mag- 
nitude limit of r* ~ 19 set by the POSS-I R plates; using 
a wider range of Galactic coordinates from future SDSS re- 
leases will help to overcome this degeneracy. 

A similar problem afflicts independent analysis of the 
density distributions of the SGS and NGS. The smaller 
ranges of Galactic longitude and distances (Figure tTTt sam- 
pled in the SGS and the decreased sensitivity to the density 
profile parameters in this direction means that neither a or 
q can be adequately constrained for the SGS taken in iso- 
lation. The NGS alone does allow limits to be placed on a 
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Figure 23. Contours of equal \ 2 constraining the spheroid den- 
sity models (a,q). Although q cannot be adequately constrained 
with these data, we find a = — 3.15± 0.30. The dashed lines show 
the one standard error interval from the best-fit a. The standard 
error in a is calculated from bootstrap resampling. 



and q, and these agree well with the values derived for the 
joint SGS/NGS sample, albeit with larger errors. 

The standard error of the a derived for the SGS and 
NGS is obtained from bootstrap resampling. The best-fit 
a is calculated for each of 1000 bootstrap samples, all of 
the same size as the original sample and chosen from it by 
randomly selecting stars with replacement. The standard 
deviation of a for these bootstrap samples is an estimate of 
the standard error in the power law index. 



7.2 The subdwarf luminosity function 

The combined luminosity functions of the SGS and NGS 
fields for the best-fit power law a = —3.15 are shown in 
Figure 1241 (the luminosity functions for our samples are in- 
sensitive to the value of q.) The results for each stripe are in 
excellent accordance with each other, with all but one mag- 
nitude bin agreeing within the la error bars. This suggests 
that there are no systematic spatial effects in the analysis of 
the two samples. 

T here is also good agr e ement with the kine matic stud- 
ies of iGizis fc Reidl Jl999t). iDahn et alJ (Il995h (scaled by 
0.75 to account forlc aseTteuio^R^!rtmu^^ 
kinematics as in lGould et al]ll99gft and lGouldl J2003l) . This 
indicates that the correction of our sample to the solar 
neighbourhood data by an r~ 3 ' 15 power law is a good ap- 
proximation, and that this distribution therefore well de- 



scribes the inner s pheroid out to 2.5k pc. The disagreement 
with the result of iGould et alJ il998l) is pronounced, how- 
ever, and this may be du e to systemat i c erro rs affecting 
their sample: for example, IGizis fc Reidl il999l) postulated 
that t he choice of local calibrating subdwarfs in lGould et all 
lll99Sf) could cause their luminosity function to be underesti- 
mated. Alternatively, there could be a real effec t behin d the 
discrepancy. The inability of the IGould et alJ |l998) data 
to match the local luminosity functions when fitted with 
a (a,q) = (—3.13,0.8) density power law indicates that 
this model does not accurately describe the spheroid. As 
we find that an r -315 power law accurately fits our data 
out to 2.5kpc, this provides more evidence that a differ- 
ence in the axial ratios of th e inner and outer spheroid 
iSommer-Larsen fc^fien|^)9C ) may be responsible for the 
discrepancy of the TGoTude^rt Il998l> result. However, there 
is growing evidence that the usual four-component model 
for the Galaxy (with discrete thin and thick discs, halo and 
bulge) is inadequate, and that a model with more continu- 
ous distributions in stellar age, kinematics and metallicity 
is mo re appro priat e JChiba fc BeenfeOOdlSiegel et al .120021 : 
lYannv et ai1l2003l) . It is possible that the deficiencies of this 
over-simplistic picture of Galactic structure is behind the 
discrepancies seen here. 

A further feature evident from Figure 1241 is the preci- 
sion of our results. That our error bars are so s mall reflects 
the size of our sample, with only the Gould (2003) lumi- 
nosity function exhibiting comparable errors . However, ou r 
subdwarf selection is more rigorous that of iGouldl i2003t) . 
who identified ~4500 subdwarfs simply by taking cuts by 
eye in the reduced proper motion plane and whose sample 
is perhaps more susceptible to thick disc contamination. 

Figure 1^21 plots the luminosity functions for stars of dif- 
ferent metallicities, with the samples divided into subdwarfs 
with metallicities either greater or less than the median of 
[m/H] x -2.4. As in the study of lGizis fc Reidl lll999l) . we 
find that the very metal-poor subdwarfs extend to fainter 
absolute magnitudes and tend to have a higher space den- 
sity than the more metal-rich ones. 



7.3 The (V/V maix ) test 

The overall completeness of the subdwarf samples can be 
estimated by using the (V/V ma x) test. For each star the ratio 
of the volume V (corresponding to its distance d) to V max 
is calculated (Equation 1301 . and the mean of this quantity 
should be 0.5 for a complete survey evenly sampling the 
survey volume. The error in this mean is 1/ (12N)? , where 
N is the number of stars in the sample. 



V 



Vmax 



d 



(30) 



The values of (V/Vmax) for each field in the SGS and 
NGS are plotted against Galactic longitude in Figure 1261 
The SGS has a combined value of (V/V max ) = 0.495± 0.010 
and the NGS has 0.506 ± 0.007. Although this is not a rig- 
orous test of completeness, especially for non-uniformly dis- 
tributed samples, these results nevertheless provide evidence 
that no significant incompleteness affects the samples. 
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Figure 24. The combined luminosity function for the SGS (solid line, open circles) and NGS (solid line, open squares) assuming the 
best-fit spher oid density power l a w a = —3.15. Th e SGS and NGS results are in excellent agreement with each other and agree w ell with 
the re sults of lDahn et alj ll995l) TCasertano et al.| dl99dfl (dot-dash, asterigks^J ^izisfc Reidl dl999l) (dashed, filled squares) and iGouldl 
J2003t) (dotted, filled triangles), but not with the outer spheroid sample of lGouldet^lT^^9a) (solid line, filled circles). This lends weight 
to the suggestion that a single power law cannot describe the density distributions of both the inner and outer spheroid. [ Note that the 
iDahn et al.lll995l luminosity function has been scaled by 0.75 to reflect use of the ICasertano eT*a l. 1990 kinematics, as in lGould et alj 
1998.] The small error bars of this st udy's luminosity function reflect the much greater size of the subdwarf sample compared to the 
other kinematic studies; IGouldl (2003) has a larger sample size, although his is more likely to suffer from thick disc contamination. 



8 SUMMARY AND FUTURE WORK 

We have demonstrated the effectiveness of the method of 
reduced proper motion at selecting spheroid stars, and have 
used it to obtain one of the largest samples of known 
spheroid subdwarfs. From this we derive the subdwarf lu- 
minosity functions to unprecedented accuracy in two dia- 
metrically opposite lines of sight in the Galaxy. 

The large samples along different lines of sight in this 
study have enabled us to constrain the form of the spheroid 
density distribution, out to heliocentric distances of 2.5 kpc, 
to a power law with an index of a = —3.15 ± 0.3. This is 



in accordance with other recent results iGould et al 


199S 


Sluis & Arnoldlll99ci lYannv et alj 


|200ri llvezic et al. 


2000 



Chen et alJl200 it ISieeel et al.ll2002l) . We are unfortunately 



unable to adequately constrain the spheroid flattening pa- 
rameter q with this study (we can only rule out a very flat 
spheroid), due to insufficient survey depth. 

Our luminosity functions agree well with other recent 



local derivations, so that with this result the solar neigh- 
bourhood subdwarf luminosity function is now well denned 
to Mv « 12.5. Our data corrected by an r -3 ' 15 spheroid 
density distribution closely match the local luminosity func- 
tions, suggesting that this power law well describes the in- 
ner spheroid density profile. On the other hand, our result 
further confirms the discrepancy between the local luminos- 
ity functions and t he outer spheroid luminosity function of 
IGould et al] dl99ct) corrected by an r -3 ' 1 po wer law. This 
provid es additional evidence that either the IGould et al] 
(1998) result is affected by systematic errors, or that the 
inner and outer spheroid follow quite different density dis- 
tributions. 

A progression of this study to larger volumes will enable 
limits to be placed on q and will provide stronger constraints 
on a. Adding the SIRTF First Look field from the EDR 
will immediately increase the sample size by ~18%, and will 
provide additional lines of sight, important for determining 
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Figure 25. The luminosity functions for stars with metallici- 
ties either greater (dashed line) or less than (solid line) the me- 
dian metallicity of \m/H\ a -2,4 in the SGS (circles) and NGS 
(squares). As in iGizis fc Reidl ll99Sft . we find that more metal- 
poor subdwarfs have a higher space density. 
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Figure 26. (V/V max ) for fields in the SGS and NGS . There is 
no evidence for incompleteness in the samples. 

a. Progressive SDSS data releases will significantly increase 
the scientific return of this work. With the first formal SDSS 
data release (DR1) some 20% of the photometry is now avail- 
able, increasing to ~50% in January 2004. This represents 
a huge increase on the 5% from the EDR, and will mean 
that tens of thousands of subdwarfs will be found with the 
methods described in this paper. These will yield even more 



accurate estimates of the subdwarf luminosity function and 
will provide much tighter constraints on the spheroid density 
law through the increase in sample size and lines of sight. 

A further development of this study that is already un- 
der way is to obtain spectra of the candidate subdwarfs. 
Confirmation of their spectral type is important for deter- 
mining the accuracy of the RPM selection methods, and 
for ascertaining the true level of thick disc contamination 
in the samples. In addition, radial velocities will enable the 
full 3D angular momentum properties of these spheroid stars 
to be investigated over large scales, with accuracy sufficient 
for the detection of spheroid kinematic sub structure (see eg. 
iHelmi et al1ll999l : iHelmi fc de Zeeuwll200d) . Spectra of sev- 
eral hundred candidates have already been obtained using 
the 6dF multi-fibre spectrograph on the UKST at A AO. 

A subsequent step to be made in the near future is to 
convert the luminosity functions into an accurate subdwarf 
mass function by use of a mass-magnitude relation. As de- 
scribed in Section this has a number of important appli- 
cations, and is particularly important to the theories of star 
formation and evolution. Future analyses will also take into 
consideration the expected binarity fraction of the samples 
in the derivation of the luminosity functions. 

This work clearly demonstrates the great potential of 
combining the old-style photographic surveys with the newer 
CCD programmes. This combination makes optimum use 
of both datasets, utilising the accurate astrometry and long 
time baselines available from surveys such as the SSS, whilst 
taking advantage of the accurate CCD photometry from 
such studies as the SDSS. 
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